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We investigate the effect of non-linear dispersion relations on the bispectrum. In particular, 
we study the case were the modified relations do not violate the WKB condition at early times, 
focusing on a particular example which is exactly solvable: the Jacobson-Corley dispersion relation 
with quartic correction with positive coefficient to the squared linear relation. We find that the 
corrections to the standard result for the bispectrum are suppressed by a factor where p c is 
the scale where the modification to the dispersion relation becomes relevant. The modification is 
mildly configuration-dependent and equilateral configurations are more suppressed with respect to 
the local ones, by a factor of one percent. There is no configuration leading to enhancements. We 
then analyze the results in the framework of particle creation using the approximate gluing method 
of Brandenberger and Martin, which relates more directly to the modeling of the trans-Planckian 
physics via modifications of the vacuum at a certain cutoff scale. We show that the gluing method 
^-», overestimates the leading order correction to the spectrum and bispectrum by one and two orders, 

respectively, in — . We discuss the various approximation and conclude that for dispersion relations 
not violating ^KB at early times the particle creation is small and does not lead to enhanced 
contributions to the bispectrum. We also show that in many cases enhancements do not occur when 
modeling the trans-Planckian physics via modifications of the vacuum at a certain cutoff scale. Most 
notably they are only of order O(l) when the Bogolyubov coefficients accounting for particle creation 
are determined by the Wronskian condition and the minimization of the uncertainty between the 
field and its conjugate momentum. 
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00 ■ I- INTRODUCTION 

CO 

The interest in non-Gaussian features of density perturbations has been mounting in recent years. It is prompted 
• , by the upcoming release of the data of unprecedentedly accurate experiments such as Planck and by more advanced 
research activities. In particular, a notable interest descends from the possible sensitivity of higher-point amplitudes 
on high-energy physics, namely the existence of new physics at high energy scales. This is due to the fact that higher 
correlation functions are suggested to be more sensitive to such high energy effects than the power spectrum itself. 

Non-Gaussianity is the non-zero value of higher-points (> 2) functions 1 of the comoving curvature perturbation 
£. In particular, here we will be interested in the bispectrum, which is obtained from the three-point function of the 
curvature perturbation. The contributions to this correlator arise from different type of diagrams depending on how 
we define the curvature perturbation (that is, if we use non-linear field redefinitions). In particular we are interested 
in the connected three-point function and in the effects that modification of the high-energy theory can impart to the 
standard slow-roll result. The outcomes of this analysis are important in two respects: for the observation of such 
effects but also for the potential risk for perturbation theory, which might break down (trans-Planckian issue). 

The three-point function for the conventional single-field slow-roll models of inflation was originally calculated in 
[l[ and found to be slow-roll suppressed. Any non-Gaussian signature is therefore a smoking-gun for deviation from 
the orthodox picture. The effect one would generally expect are two-fold: 

• a peculiar shape function 2 

• an enhancement factor for some specific configurations. 
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1 The two-point one is related to the power spectrum. 

2 That is a peculiar dependence of the result on the external momenta, leading to a particular shape for the graph of this function. 
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Generally these effects are combined, and we have particular enhancement factors when the external momenta of 
the three point function (in momentum space) have specific configurations. 

Various methods have been employed to model high energy modifications to the standard theory. For example, the 
effects of the trans-Planckian physics have been implemented via a stringy-modified space-time uncertainty relation 
0, or utilizing the concept of a boundary action, which provides the boundary conditions for the perturbations fields, 
subject to the renormalization from higher energy scales 0, 0] ■ A particularly interesting method is related to the 
choice of boundary conditions (the so-called "choice of vacuum" ) for the solution of the equation of motion of the 
Mukhanov variable at the new physics hyper-surface (NPHS), which corresponds to the time the physical momentum 
of perturbations reaches the scale A i$r7\- In this case, the power spectrum exhibits superimposed oscillations whose 
amplitude is given by H/A, where H is the Hubble parameter during the inflation. The bispectrum is instead found 
to be modified by enhancement factors for the enfolded/flattened configuration of momenta, whose real magnitude 
depends however on the value of the Bogolyubov parameters calculated within this approach (8l-[l0|. In such a 
configuration, two of the momenta are collinear with the third one in the momenta triangle. 

In an alternative approach to model the effects of trans-Planckian physics, that we will mainly focus on here, the 
equation of motion for the perturbations are modified and new dispersion relations, which differ from the linear one for 
physical momenta larger than a fixed scale of new physics p c , are considered (T^ |. Such modified dispersion relations 
represent the violation of Lorentz invariance and were also employed to analyze the possible effect of trans-Planckian 
physics on black hole radiation [l3j. Modified dispersion relations of the kind we will investigate here, can also be 
derived naturally from the recently proposed Hof ava- Lifshitz gravity dUU] which is a non-relativistic renormalizable 
modification of the gravity in the UV region. They also arise in effective theory of single field inflation when the 
scalar perturbations propagate with a small sound speed [l7[ . 

In this paper we focus on examples where one can obtain the exact solutions to the field equations and therefore has 
control over subtle effects. As the evolution of the mode function after the modified dispersion relation has entered the 
linear regime can be mapped to an excited state in the NPHS approach, using the exact solutions we have a precise 
understanding and quantification of particle creation and can understand the reason that lies behind the presence or 
absence of the enhancement factors in the bispectrum. The general result that we find is that, in absence of violation 
of the WKB condition at early times, the enhancement factors are not present in the bispectrum. On the other hand 
one would expect a more favorable situation for dispersion relations with early time WKB violation. However, as we 
have found no example of such modified dispersion relations where the field equation was exactly solvable, we discuss 
this case only briefly and leave the analysis of this point to future research. 

We also show how the enhancement factors are at maximum of order O(l) in some of the emergent approaches to 
trans-Planckian physics in the NPHS framework. Notably, this occurs in the case where the Bogolyubov coefficients 
accounting for particle creation are determined by the minimization of the uncertainty relation between the field and 
its conjugate momentum and by the Wronskian condition [5], which yields the largest correction to the spectrum. 

We will first study a minimal modification to the linear dispersion relation, known as Jacobson-Corley (JC) [l3|, 
both solving the field equations exactly, and via approximation methods. The JC dispersion relation has a positive 
quartic correction to the linear term and there is no WKB violation at early times 3 . 

We find no modulation in the modified spectrum of perturbations, and the amplitude of perturbations is damped 
quadratically with the increase of the ratio — . This is explained in terms of the very small particle creation, due to 

Pc 

the absence of WKB violation at early times. 

We then turn to the analysis of the modification to the bispectrum: using the exact solution to compute the three- 
point function, we find no large enhancement factors, and no interference terms as the one studied in |8l4l0|. We then 
try to understand and inter pret these results in the language of particle creation using the gluing approximation of 
Brandenberger and Martin [12[ for the solutions. We show that the gluing method overpredicts the magnitude of 
modification. 

The outline of the paper is as follows: we review the formalism, the notation and the techniques to solve the field 
equations in section |Hj We then study the JC modified dispersion relation in section Mil first we solve exactly the 
field equations in section IIII Al then we study the two- and three-point functions in sections IIIIBI and IIII CI We 
then try to explain and interpret the results in the framework of particle creation, and make contact with the NPHS 
framework in section IIVI We finally discuss the general results, present the outlooks regarding modified dispersion 
relations with early time violation of WKB condition, and conclude in section FVl 



3 The only WKB violation occurs when the modes exits the horizon at late times, when the dispersion relation has become approximately 
the standard linear one. 
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II. FORMALISM AND NOTATION. 



The starting point for our analysis is the definition of the curvature perturbation of the comoving hypersurface 

d 3 k 
(2^0 

where Cviv) is defined in terms of the mode function uArf) as 



CM = / 7^TC^)e ifc ", (1) 



CM = (2) 



with equation of motion 



uZ+(u(r ] ,k) 2 --)ut = (3) 



z 



where 4 z = ^ and oj(r),k) is the comoving frequency as read from the effective action. In the standard formalism 
uj(r],k) is taken to be equal to k = |fe|. Instead, a modified dispersion relation would correspond to considering 
different dependence of ui on fc, r\. For isotropic backgrounds, which we will confine ourselves to, one expects uz(r\) 
and ui{rj, k) to depend only on fc. Hence, we may drop the arrow symbol at some points. 
The three-point function is given by the formula 

(C(zi)C0caX(s3)) = -2Re (J" d^(^n|C(^)C(^)C(^3)^i(r?')lV'm)J (4) 

where Hj = — J d 3 xa 3 (-^) 4 H( 2 d~ 2 (' is the interaction Hamiltonian as defined in ri- ln and \ipi n ) are the initial 
time and state (vacuum). The standard result in slow-roll inflation is obtained with rj- m = — oo, and \ipin) is taken to 
be the Bunch-Davies vacuum. 

In momentum phase space the three-point function becomes [l[ 

(CkM<k 2 (vK%M) = 2Re (- l ( 2 ^ 3,5 (E^)(|) 4 ^ f n drf^-fid^M) + permutations^) (5) 

where rj is taken to be a late time when all three functions are outside the horizon 5 . (77, 77') is the Whightman 
function defined as: 



^ _ H 2 u fc fa) u%(rf) 
<j> 2 a(v) a(v') 

We define the bispectrum following [l8| 



(^MC^vXnM =S(J2^)(2n) 3 F(k 1 ,k 2 ,k 3 ,r ] ), (7) 

where the translational invariance has imposed the conservation of momentum. Scale-invariance requires that the 
function, F, to be a homogeneous function of degree —6 and the rotational invariance imposes it to be only a function 
of two variables, say x-i = kijk\ and £3 = ks/k\. To avoid counting the same configuration twice, it is further assumed 
that X3 < X2- The inequality 1 — X2 < £3 also comes from the triangle inequality. 

The function F contains a lot of information about the source of non-Gaussianity and could be used to distinguish 
among different models. The limit in which the configuration is such that X3 ~ and xi — 1 is recognized in the 
literature as the local configuration. The one in which x<i ~ X3 ~ 1 is known as equilateral one. 

Inflationary models in which non-linearity is developed beyond the horizon, tend to produce a more local type 
of non-Gaussianity. On the other hand, for the ones that the correlation is among the modes with comparable 



wavelength, equilateral type of non-Gaussianities tend to be produced 18jj. These modes will exit the horizon around 
the same time. 



4 4> being the infiaton field. 

5 In the case of modified dispersion relations, at this time the relation has become effectively the standard linear one. 



4 



A. Solving the field equation. 

In the case of a modified dispersion relation, equation ([3]) is often solved by approximation [l2[ . Two methods have 
been used: one is the WKB approximation. This approach is convenient, as it is independent of small parameters in 
the dispersion relation, because it only requires a slow variation of U 2 = lo 2 — — in time, i.e. [20]. 



Q 




U" 


u 2 




2C/ 3 



3^ 



« 1 (8) 



The other is the gluing procedure of Brandenberger and Martin, where approximated solutions valid in different 



regions of behavior of u) are glued together asking for the continuity of the functions and their first derivatives [12j . 

In the NPHS approach, enhancement terms appear in the bispectrum due to the interference terms between negative 
and positive frequency parts of the Wightman function as one starts from a non-Bunch-Davies vacuum. In the case 
of modified dispersion relations, similar interference term will generally arise solving Q using the gluing procedure. 
However, in all cases the actual presence/absence of enhancement depends on the magnitude of the second Bogolyubov 
parameter /3g [8l4l0f. 

With the knowledge of the exact solution in some specific examples, one can gain more precise information about 
the shape of the bispectrum. We want to explore these features in the case of modified dispersion relations and in 
particular verify the presence/absence of enhancement factors and their difference with the NPHS case where one 
employs the cut-off. To do this, in the following we will consider a particular (but very illustrative) example where 
the exact solution can be found. We then analyze and interpret the results in the framework of particle creation with 
the approximation methods described above, and discuss the differences. 

III. CORLEY-JACOBSON DISPERSION RELATION WITH bi > 0. NO VIOLATION OF WKB 

APPROXIMATION. 

In this section, we will focus on Jacobson-Corley dispersion relation with positive quartic correction. We assume 
that at very high physical momenta, p > p c , the linear dispersion relation lo — p gets modified as follows: 

u 2 phys = F 2 (p)=p 2 + b -^- (9) 

Pc 

The dispersion relation is motivated by studies in black-hole physics (l3j |. and was later used to study the effect 
of trans-Planckian physics on cosmological perturbations [lU, [l9j]. Here, we will find an exact analytic form of the 
solution and will derive the effect of this modified dispersion relation on the two and three point functions analytically. 
Then, we study the same problem making use of approximation techniques and discuss the results. 

A. Exact solution of the field equations 

In usual quantum field theory, the cosmological perturbations satisfy eq. (|3J) as the equation of motion. A modified 
dispersion relation of the kind of Corley-Jacobson entails the replacement 

k 2 =► a 2 ( V )F 2 ( * ) = k 2 + b,^- (10) 

Then in de-Sitter space-time, where a = —77^, equation ([3]) reads: 

<(r,) + (fc 2 +AV-iH = (ff) 
V 



where, 



^. (12) 
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The above differential equation has two independent closed-form exact solutions that can be given in terms of 
WhittakerW(a, b, z) and its conjugate WhittakerW* (a, b, z) functions. For brevity, we will show these two functions 
as WW(a, b, z) and WW* (a, b, z) from now on, respectively. The exact solution is then: 

Ci I 1 3 • , 2 2 \ , c 2 ww * / i 3 2 2 



M fc = —WW — , -, -kfcV + ^=WW* — , -, -ieifcV 

Cl tt™V « 3 ■ ,2 2 \ , C 2 3 , , 



WW — , -, -icfcV + -7=WW — , -, iefcV . (13) 



The above solution is subject to the Wronskian condition 

u{rf)u'*{rf) - u*(rf)u'(rf) = i, (14) 
which impose the following constraint on C\ and C2: 

2i|Ci| 2 exp( — )ek 2 -2i\C 2 \ 2 exp( — )ek 2 =i, (15) 
4e 4e 

This condition would not uniquely determine C± and C 2 . One has to make extra assumptions to determine C\ and 
C-2- We will assume that the mode function approaches the positive frequency WKB at ealy times. As it was shown 
by [Hj], this choice will minimize the energy density too. 

In the limit of rj — > —00, the exact equation of motion is reduced to the following form 

u'Hr,) + e 2 kW<r,) = 0, (16) 

and its approximate positive frequency WKB solution is: 

«*(„) * -jL=^{-i f U (r,)drf) = ^^exp(^). (17) 

The limit of the exact solution (| 13[) for large time can be found writing the Whittaker functions in terms of KummerU 
functions and using the asymptotic form of the latter for large argument [2l[ . We obtain 6 

uM = -^=ex P (^)H e fcV) 5/4 u(^,^,-^v) 

Cl = (-iefcV)* exp(^-^-) + -^L(iefcV) - * exp(-^-^-), 



(18) 

where we see that the positive frequency WKB solution corresponds to the choice C2 = 0. Using the Wronskian 
condition, (fTo"]) . G\ is fixed up to a phase: 

C X = (19) 
V2e/c 



so that our exact complete solution to the mode equation is, finally, 



The extra factor of (— iek 2 rp) ^ is the subleading correction to the WKB approximation in the limit of rj — > —00 
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FIG. 1: The figure shows the dependence of the power spectrum in de-Sitter space to the parameter e, which quantifies the 
ratio of Hubble parameter over the momentum scale p c around which the behavior of the dispersion relation changes. 



B. Two-point function and power spectrum 

With the exact mode function in equation (fT5|) one calculates the power-spectrum 



Pn{k,e) 



k^lP 
2^1^ 



u k (v) 



a(rj) 



(21) 



The spectrum is found to be still scale-invariant (no dependence on k), but the standard result now shows corrections 
depending on e: 



Pn(e) 



exp(-£) 



9 16^/ 2r( | _ £ )r( | + £ } 



(22) 



The plot of the power spectrum in terms of e is shown in figure (JXJ) . The scale- invariance of the power spectrum is 
again due to the scale-invariance of the de-Sitter background. The corrections to the standard slow-roll are 



These corrections are appearing at second order in e. 



H 2 

2nd) 



5 2 

1 - -e 2 

4 



(23) 



C. Three-point Function 

Given the exact mode- function (f20|) . the Wightman function for the positive frequency WKB vacuum is 



2ek 2 



i 3 
4? 4' 



iek 2 r] 2 ] WW 



-i 3 
47' 4' 



2„'2 



iek i] 



(24) 



To compute the three-point function from equation (|5|) , the first argument of the Wightman function has to be set at 
a time when the mode k is outside the horizon, generically taken to be time zero. Then, the Whightman functions 
must be differentiated with respect to the second argument. We obtain 



a(r,)d v G>(0,v)) = - 



(-l)itf 3 exp(-£)^ 



e ffci^(-„)*r[f-£] 

x WW ( — ,-,iek 2 n 2 
V 4e 4 



4< WW ( 1 - JL, |, iek 2 r] 2 ) + (i + e + 2ie 2 k 2 V 2 ) 



(25) 
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For computing the three-point function, it is convenient to write the WhittakerW functions in terms of KummerM's 
based on the following relation which is valid in the principal branch, — tt < arg(z) < tt. 



WW(a,M) = exp(--) 



/ r(26)M(i - a - b, 1 - 2b, z)z^- h r(-26)M(± - a + b, 1 + 26, z)z^ 



r(i-a + 6) r(i-a-6) ) 

We can then exploit the following expansion of the KummerM functions in terms of Bessel- J functions (2lj : 



(26) 



M(a, b, z) = T(b) exp(hz) V C n z n (-az)^ 1 - b - n ^ J b _ 1+n (2^f^z) . 



71=0 



C„+i is given in terms C„, C„_i and C„_2 from the following recursive relation 
1 



n+l 



(n + 1) 



[((1 - 2h)n-bh)C n + ((1 - 2h)a- h{h - l)(b + n- l))C„-i - h(h - l)aC n _ s 



Co 



Ci 



-bh. 



c-> 



l -(2h-l)a+ h -{b + l)h 2 



(27) 



(28) 



(29) 



Although strictly valid for small z, the above expansion becomes the correct expansion of the bispectrum for small e 
when the integration over time in the formula ([5]) for the three-point function is performed. We will also prove this 
point in section HV 11 

As it appears from equations (|25|) and (|26|) . the function a(r])d v G^ (0, ij)) is proportional to exp( ~ ze ^ v ). Factoring 
out this term, one can expand the coefficients in powers of e after the KummerM functions are substituted with their 
approximate Bessel series expansion. To second order in e, we have: 



H 



-iek Tf 



2„2 



ikrj 



1 

2k 



ikrj 2 



5 5t7 I % I 

i— kn 2 k 2 n 3 H fc 3 ?7 4 

8fc 8 8 ' 12 ' 16 



(30) 



It is now straightforward to obtain the three-point function. We will perform the relevant integration over 77 and 
expand the final results in series of e in order to obtain the corrections to the standard (e = 0) spectrum. We show 
this procedure in detail for the zeroth order term in the square bracket of (I30[) . we then present the final results 
regarding the other terms in the square brackets. 

The zeroth order term in the square bracket of (|30[) leads to the contribution 



3x3 / 




Mp 2 H 



1 



fe 3 



c.c. + permutations, 



with solution 



(CjsjCjgCjg)o 



b e k 1 k 2 k 



■ exp 



-ir)(ki 



H° 



(-1)271-2 exp(-^)Erfi 



(l-i)(-fct+efcjt)) 



b 2 M 2 



Ti=0 



permutations, 



(31) 



(32) 



where 3? indicates the real part and 



k t = h + k 2 + h , k 2 = kf + k 2 + kj 



(33) 



In contrast to the Lorentzian dispersion relation, u> — p, the integrand remains finite at r\ = —00 even without the 
change of variable from r\ — > r\ + «e| 77] . In this limit the Erfi function tends to i. This is interesting, as the adhoc 
prescription to make sense of the contribution at —00 is completely taken care of by the modified dispersion relation, 
which is expected from quantum gravity. 

Taking the real part and expanding to second order in e, one obtains 




2H TT 



6^3 



PldkikaktM* 



6H 6 k 4 s TT 3 



1 



1 

k 2 



1 



(34) 



The first term in the square brackets is the regular quantum field theory result in absence of any nonlinearity in the 
dispersion relation. The first correction is proportional to e 2 . 

We then apply the same procedure of integration and expansion for small e for the higher order terms in the square 
brackets of (|3U|) . We obtain 
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AF(1,x 2 ,x 3 ) 
F(1 , x 2 , x 3 ) 




x 2 



FIG. 2: The figure shows the dependence of ^S 1 '* 2 '* 3 ? on xi and 3:3 setting e = 0.1 



for the term of order e in the square brackets 



- for the term of order e 2 in the square brackets 
where 

fee = fc? + *a + fc f ■ ( 37 ) 
Having all the contributions, we can compute the total three-point function: 

(fefefe)tot = (fefefe)o + (fefefe)i + (fefefe>2, (38) 



AF(fct^,"fc 3 ) 

relation, is 



and the relative change in the shape function F, L> _j 3 to second order in e, due to Jacobson-Corely dispersion 

F(k 1 ,k 2 ,k 3 ) 



AF(k[,t 2 ,t 3 ) _, . , fee 3 fc 2 ^ 2 
F{ki,k 2 ,k^) 



M+t|-^K. (39) 



We have plotted AF (^i^°- 002 ) m hV ur e[21 setting e = 0.1. As one can see from the plot and the formula (1591 . the 

F(fei,fe 2 , 0.002) t> i_i t> f " — i" 

general effect up to quartic order in the expansion is very mild and proportional to e 2 . The modification is slightly 
configuration-dependent and the equilateral configurations are more suppressed with respect to the local ones by a 
factor of one percent. One should also notice that the enfolded configurations are not enhanced with respect to the 
other ones. 



IV. EXPLANATION OF THE ABSENCE OF ENHANCEMENT FACTORS IN THE FRAMEWORK OF 

PARTICLE CREATION 

We have seen in the previous section that no interference term seems to appear in the bispectrum. This seems to 
contradict the expectation from the fact that a modified dispersion relation can be understood in terms of a "vacuum" 
state for the perturbations given by an excited state. 

Comprehending this point is crucial in showing that there will not be any enhancement for the enfolded config- 
urations. In particular, the analysis arriving at this result is based on an expansion of the Whittacker function in 
(|27p valid for small z = iek 2 r/ 2 : as the integral in the three-point function ([5]) goes up to 77 — > —00, we would like to 
understand in other terms the physical reason behind the absence of the interference terms. 

To make sure that it does not depend on the particular expansions used, we have performed another analysis 
without using any series expansion. Finally, we have also considered the approximated solution found with i) the 
gluing method, ii) the WKB approximation and compared the different approaches. 
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1. Analysis with complete solution in the region u) p h ys > H 



In this section we are going to solve the equation exactly in the region w p hy S > H. The solutions are more 
manageable than before and need not be expanded. We make a change of variables: x = —V^ekrj, so that (fTTj) 
becomes 



d l^ + ^ + Y e -^ Uk = Q - (40) 

We see that for \ 2 > 4e O k\q\ > 2, we can approximate the equation as 

d 2 x u k + (^ + ^)u k =0. (41) 
The general solution to this differential equation is 

ufe(x) =c 1 E(a,x) + c 2 E*(a,x), a =~^' ( 42 ) 
where E{y,y) are complex parabolic cylinder functions (2~i| . 

2 2 

The solution has two asymptotic regions, A and B, corresponding to < \a\ and ^- > \a\. Using the property of 
parabolic cylinder functions |2l|, one finds in A 

E(a, X ) ~ , 1; 2 1 4 , -y (V2ekr,)i e i *'?#+ i * +i *> , (43) 

X 2 »M 64 fc2 |7j| 2 

where (f>2 = arg(r(i — j-)). We reproduce the correct WKB asymptotic (see [l9|) choosing 

e — if — i<t>2 

ci = ' 5 x, j , c 2 = 0. (44) 
In region B we use the asymptotic of S(a, x) for X ^ °j obtaining 

e *f e -ik n +0(e\, 2 ) ( 45 ) 



1 

2 4 



r(!-^) 



We therefore find that no term proportional to e lkv appears in the mode function, itfc. Since the Wightman function 
is proportional to u' k * , we are assured that no interference term will appear in the Wightman function. We will now 
make a comparison with the result obtained using the gluing and the WKB approximations. 



2. Analysis using gluing and WKB methods 

In this section we approximate the solution using the gluing method. We also discuss the WKB approach. As for 
the gluing, to start with we notice that there are three different asymptotic regions for equation (jlll) and its solution: 

(fc|?7|) 2 < 2 region I 

2 < (k\ij\) 2 < (T 2 region II 

{k\r)\) 2 > c~ 2 region III 

(46) 



One approximates the solution in the various regions in (|46|) and asks for the continuity of the mode functions and 
their first derivatives across the boundaries of the regions. We start from region III, where, being k\rj\ > e _1 , equation 
(1111) can be approximated by (|16p . The asymptotic solution in region III is given by (|17p. where we have taken only 
the positive frequency mode. 

In region II we have instead a solution of the form 

u k = 4Le- lfc " + iLe 4fc " . (47) 
V2k V2k 
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The coefficients a k , (3k are fixed imposing the continuity conditions at the time rj c given by 



which separates regions II and III. We find 



* = -i> (48) 



(49) 




where we have used k\rj c \ — e _1 . We can now investigate the appearance of enhancements, expected from the 
interference terms that would appear in ([5]) due to the form of the solution (|47|) . In fact, the integrand of the 
bispectrum given by ([5]) is proportional to the product of three Wightman functions. The integral domain over rj 
can be divided into regions I, II and III where we can use the approximate solutions in each region. In particular, 
because of the form of the solution in region II, see (|47p. in parallel with the results of [9j one finds enhancement 
factors multiplying the rest of the bispectrum from the interference terms. However, from (|49|) these factors are of 
the order of 



\(3 k kn c \ ~ e\k Vc \ = 1. (50) 

That is, there is no large enhancement factor, even before taking into account the suppression that comes from the 
2D projection on the CMB surface. 

Furthermore, this method also overestimates the corrections to the power spectrum. This can be seen by comparing 
to the WKB-method, where there is no distinction between regions II and III. By matching the WKB solution, valid 
there, with the growing and decaying modes solving the equation of motion in region I, one readily finds that there 
are no interference terms, and that the spectrum is corrected only at the order e 2 , in agreement with the result 
obtained with the exact solution. Failure of gluing method in obtaining the correct amount of modification to the 
power spectrum was encountered in [22] too. 



a. Contact with the NPS method and importance of the Wronskian condition 

The NPHS procedure for modifying the vacuum in presence of trans-Planckian physics makes direct contact with 
the result we have found in this section via the gluing method, if we identify the cutoff with the scale p c . In that case, 
the evolution in region III is accounted for by an excited state implementing boundary conditions for the perturbations 
that lead to the conditions (|4"9")l at the k— dependent junction time (|4"8)> . 

Also in the NPHS scenario where the Bogolyubov parameter is determined by the minimization of the uncertainty 
relation between the field and its conjugate momentum and by the Wrosnkian condition, there is actually no large 
enhancement factor. This descends from the fact that in that case [5J 

\Pk\(NPHS) = (SNcI)- 1 ~ ^ = e (51) 

where we have used the definition (|48[) . Furthermore, in other prescriptions for the choice of vacuum in line with the 
NPHS approach [jjl], the coefficient @k is even further suppressed, as powers of Therefore, it seems that in most 
of the emergent cases in the cutoff procedure the second Bogolyubov factor is very constrained and will not lead to 
any large enhancement factor. 



V. DISCUSSION, OUTLOOKS AND CONCLUSION 



Quantum gravity effects are expected to modify the standard Lorentzian dispersion relation at very high energies, 
close to the Planck scale. One notable example is Hofava-Lifshitz gravity [lj, EH, where higher order derivative 
corrections yield a modified dispersion relation. Such high energy modifications could leave detectable signatures in 
the temperature fluctuations of the CMBR, allowing an experimental detection of this aspect of high energy theories. 

Using the illustrative example of the Corley-Jacobson (JC) modified dispersion relation with positive coefficient, 
we were able to investigate the modification to the bispectrum in the case where there is no violation of the WKB 
condition at early times, and to compare these with the expectations deriving from the NPHS approach. The field 
equations for the JC dispersion relation could be solved exactly, so that the complete formula for the bispectrum, in 
series of the ratio e = can be obtained. In this ratio, H is the Hubble parameter during inflation, while p c is the 
scale at which the modification to the linear dispersion relation become important. 
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It has been found that the leading correction to the standard bispectrum is suppressed by a factor e 2 and no 
configuration of momenta leads to compensating factors that enhance the correction. The modification is slightly more 
pronounced for the equilateral configurations by a factor of one percent with respect to the local ones. The analysis of 
the bispectrum in two different approximations and in the framework of particle creation leads to the conclusion that 
if the modified dispersion relation does not violate the WKB condition at early times, the particle production is too 
small to generate large modifications to the bispectrum. In particular, using the WKB approximation gives results in 
qualitative and quantitative agreement with those obtained with the exact solution, while the gluing method proposed 
by Brandenberger and Martin, which relates more directly with the NPHS approach to trans-Planckian physics, is in 
partial disagreement. 

One would expect more favorable possibilities when the modified dispersion relations do violate the WKB condition 
at early times. However, this point needs to be verified, as the real presence of enhancements could depend strictly 
on the form of the dispersion relations. We leave this to future investigation. 

We have also shown that in most cases within the New Physics Hyper-surface approach to modeling trans-Planckian 
physics, where at the scale p c of new physics the theory is cutoffed and a boundary condition is imposed on the fields, 
there is no occurrence of large enhancement factors in the bispectrum, not even for the enfolded configurations studied 
in Q. Most notably, the enhancement is only of order O(l) when the Bogolyubov coefficients accounting for particle 
creation are determined by minimizing the uncertainty relation and by the Wronskian condition, which yields the 
largest correction to the power spectrum. 
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